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We present a theory for the spin-Peierls transition in 
CuGe03. We map the elementary excitations of the dimer- 
ized chain (solitons) on an effective Ising model. Inter-chain 
coupling (or phonons) then introduces a linear binding po- 
tential between a pair of soliton and anti-soliton, leading to 
a finite transition temperature. We evaluate, as a function of 
temperature, the order parameter, the singlet-triplet gap, the 
specific heat, and the susceptibility and compare with exper- 
imental data on CuGeOs. We find that CuGe03 is close to 
a first-order phase transition. We point out, that the famous 
scaling law ~ <5 2//3 of the triplet gap is a simple consequence 
of the linear binding potential between pairs of solitons and 
anti-solitons in dimerized spin chains. 



I. INTRODUCTION 

With the discovery of the first inorganic spin-Peierls 
compound CuGeOaEJ it has become possible to investi- 
gate the physics of the spin-Peierls transition in quasi 
one-dimensional spin chains with high precision. Of par- 
ticular interest has been the study ofJthfi magnetic excita- 
tion spectrum by neutron scatteringBel and Raman scat- 
tering experiments™!! It has been observedoQ, that the 
spin-Peierls gap, i.e. the gap to triplet excitations out of 
the singlet ground-state, has a temperature dependence 
which is difficult to explain with existing theories of the 
spin-Peierls trajisitiorjQ'EL 

Existing theories of the spin-Peierls transition of spin- 
1/2 Heisenberg chains are based on the mapping of the 
spin chain to an interacting gas of spin-less fcrmions via 
the Jordan- Wigner transformational! The chains are 
coupled to a three-dimensional phonon mode and the 
transition occurs when the respective phonon frequency 
becomes soft. In the spin-Peierls state the unit-cell dou- 
bles due to the alternating displacements of the magnetic 
ions. The effective Hamiltonian for the magnetic excita- 
tions along the chains in the dimerized state is thcr£3 

H = JJ2 [(l + ^-l^jSi-Si+i+aSi-Si+a]. (1) 

i 

The nearest neighbor (NN) exchange has alternating 
strength J(l ± S) due to the doubling of the unit 
cell below the spin-Peierls transition temperature T$p. 
Also included in ([!]) is a frustrating next-nearest-neighbor 
(NNN) term ~ J a. The ground-state phase diagram of 
(|lj) has been mapped, out by a density matrix renormal- 
ization group studyO. There is a gap everywhere in the 



phase-diagram except on the line 5 = and a < a c , 
with a c « 0.2411. A vividly debated question of key in- 
terest is therefore the actual magnitude of the strength of 
the frustration a. For CuGeC>3 the parameters hawe-beeti 
estimated as J w 150-160K and a » 0.24-0.35E3By. 

A consequence of existing theories of the spin-Peierls 
transition is the occurrence of a soft phonon mode abtpe 
Tsp, which has not yet been observed for CuGeOsli-X 
It is therefore of interest to pursue the possibility of a 
purely electronically driven spin-Peierls transition, as it 
could occur for a > a c . In such a scenario the phonons 
would just follow the pre-formed electronic dimerization, 
i.e. the formation of NN singlet pairs. Here we pro- 
pose a simple molecular-field type theory for the spin- 
Peierls transition by mapping the low-energy solitonic 
excitations of dimerized spin chains on an effective Ising 
model. We emphasize that this procedure is valid both 
for the phonon-driven and for the electronically-driven 
spin-Peierls transition. We present results of the mean- 
field theory in comparison with experiments on CuGe03. 
We find that a straightforward determination of the frus- 
tration parameter a is not conclusive and that CuGe03 
is close to a first-order phase transition. We also-pote 
that the famous ~ <5 2 / 3 scaling of the triplet gapae is 
a simple consequence of the linear binding potential be- 
tween solitons and anti-solitons in dimerized spin chains. 



II. MAPPING ON AN EFFECTIVE ISING 
MODEL 

A strictly one-dimensional system shows no long-range 
order at finite temperatures due to the solitonic excita- 
tions. The appropriate solitons in a spin-Peierls state 
are domain-walls in between two different dimer cover- 
ings, see Fig. |l|. In this picture a dimer consists of a 
NN pair of spins in singlet state. A single soliton in an 
otherwise (dimer-) ordered chain is inpteality a compli- 
cated objecttj. It is spatially extendedE^I and has a spin 
degree of freedom together with a disoersion, which in 
the Majumdar-Gosh model ( a = L/2 ED) has the form 
J (5/4 — cos2/c). It has been shownta, that the solitons 
yield an accurate description of the magnetic susceptibil- 
ity and hence of the Hilbert space. 

The dispersion of the solitons might be approximated, 
in general, by 



E(k) = y/(E s ) 2 + (cJ sinfc) 2 



(2) 



where E s is the gap to solitonic excitations, J the 
exchange integral, c a constant of order of unity, and 
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k is the wavenumber in units of inverse lattice spacings. 
Only the low-energy solitons are effective in destroying 
the long-range dimer order, and we approximate here (g) 
by a constant E(k) = E s . We will furthermore not 
discuss the consequence of the spatial extent of the soli- 
tons and we will include their spin-degree of freedom only 
when determining the susceptibility. We take one of the 
two possible dimer configurations of the linear chain as 
the reference state and label every "good" dimer by +1 
and every bond between two "wrong" dimers by — 1; A 
bond between a soliton and a wrong dimer is also labeled 
— 1 as shown in Fig. |l|. Every possible soliton configu- 
ration is therefore mapped to a Ising-spin configuration 
living on every second bond. A similar mapping has been 
used recently by Mostovoy and Khomskiita. 
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FIG. 1. Phenomenological picture of copper chains in 
CuGeOs,. a): A soliton and an anti-soliton in a dimerized 
chain with the corresponding values of the dimer operators 
an. There is one dimer operator for every pair of sites, b) and 
c): Illustration of two different dimer configurations which 
lead to an inter-chain coupling contribution to the energy of 
±0.5A 9 . The origin of the coupling may be attributed to the 
inter-chain exchange as well as to phononic effects. 

We now consider a two-dimensional array of chains and 
define by er^ — ±1 the Ising variable on the chain, 
where the site index i runs over every second bond of 



in the strongly anisotropic limitta. We may therefore 
decouple the chains via 



0~i,l (o~i,l+l) +CT-i,i+l {(Tij) - (ctu) (cu+i) 



(4) 



where the {ci.i) denotes the thermodynamic expecta- 
tion value. Translational invariance perpendicular to the 
chains yields (<7i,z+i) = (cy) = (a) and the Hamilto- 
nian for a single chain (the number of Ising variables is 
L/2) becomes in this mean- field approximation 

1 L/2 L/2 L L 



i=i 



(5) 



where we have set B — A g (cr). This is just the Ising 
Hamiltonian for a ferromagnetic spin chain in an external 
magnetic field B. Alternatively we can interpret (||) as 
an effective model for domain walls in a spin-Pcicrls state 
with a linear binding potential V(x) = B ■ (x + 1) (with 
x > 1 being the distance in sites, not bonds) between 
a soliton and an anti-soliton. 

So far we have neglected the spin-phonon coupling X sp . 
A finite value for the spin-singlet order parameter (a) 
leads through X sp to a finite lattice dimerization and 
consequently to a modulation of the exchange constant, 
J(l ± 5), in ([j]), with S ~ X sp (a). The spin-phonon 
coupling therefore adds a term ~ X sp (a) to the confining 
potential B. On a mean-field level this corresponds to a 
rescaling of the coupling constant 



A g -> A g + ciX s 



B 



B + c 2 S, 



(6) 



with appropriate constants c\ and ci- 

The partition function Z can be obtained from (j^) 
by the transfer matrix method. The free energy F = 
— fc^Tln Z is given in the thermodynamic limit (L — > oo) 
by 



F 

L/2 



-k B T\n\ a + -A CJ {af 



where 



a chain of length L. Including a coupling A g between Aq = C0S h(/3A g (o-)) + Jsmh 2 {f3A g (<j)) + c 
dimers on NN chains we have v 



(7) 



(8) 



H. 



2D 



E* X] 2^ _ a i,l a i+l,l) 
i,l 



1 



with (3 = l/(ksT) being the inverse temperature. From 
(3) the free energy all physical quantities can be derived. 



as the effective Ising Hamiltonian. The inter-chain cou- 



by the inter-chain ex- 
In this case A„ ~ 



pling A 9 might be induced 
change coupling, Jl ~ 0.1J. 
(J±) 2 /J. An indirect coupling of the chains via phonons 
is also conceivablea. It has been shown recently, that a 
molecular-field decoupling of inter-chain interactions in 
quasi one-dimensional systems is a good approximation 



III. EXCITATION ENERGIES 

We have two relevant excitation energies, the gap to 
single soliton excitations, given by E s in (^|), and the 
gap to triplet excitations, A^v, as measured in a neutron 
scattering experimentu. Let us now discuss the relation 
of Ajy to E s . 
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A triplet excitation can dissolve into a soliton/anti- 
soliton pair. The linear binding energy V(x) = B-(x + l) 
in between a soliton and an anti-soliton ip-i(||) leads to 
a confinement of soliton/anti-soliton pairsEI In (||) we 
have neglected the kinetic energy of the single solitons 
(||), which is of order J. The energy levels ofr-a,rparticle 
in a linear confining potential are well knownMEa. The 
lowest eigenstate has, in the limit B« J, the energy 



2E„ 



•'"7 



2/3 



(9) 



with d w 2.33. Eq. (|9|) gives then the gap to triplet 
excitations as a function of soliton energy E s and the 
strength of the confining potential B. The mean exten- 
sion of a soliton/anti-soliton pair scales like (B / J) -1 / 3 . 

For a < a c we have B — > C28 and Eq. (||) takes then 
the form 



A 



N 



2E S + 2.33 v s 



a<a c 



C 2 8 
V s 



2/3 



(10) 



,is the a-dependent spin- 
: that the gap scales for 
for 



where v s « (jff/2)(l - 1.12a) 
wave velocity^ 3 ! It is knownc 

a < a c like ~ S 2 ^ 3 implying that E s 8~ 2 /2-. — > 
5 — > 0. Comparison with numerical resultstil leads to 
C2 « 0.85. The functional form of the dependence of the 
soliton excitation energy E s on 5, or alternatively on 
the dimerization order parameter (a) is not known at 
present. It might be extracted for a < a c , in princi- 
ple, from a sub-leading scaling analysis of the excitation 
energy 



A N = dJ(S) 2/3 + 2E S (S) 



(11) 



but this has not yet been done. We have therefore de- 
cided to assume the functional form 



E s {(a)) = £oo + (E - £oo)(a} 2 



(12) 



where Eq is the zero-temperature value of E s and where 
Eoo is the soliton energy in the disordered phase, i.e. for 
T > Tsp- Eq. (12) is, in the spirit of a Landau functional 
(|l3|), the simplest form consistent with the symmetry of 
(a). 



IV. LANDAU EXPANSION 

We may expand the free energy (|^) in powers of (er), 

F = F +a(T)(a) 2 +6(T)(a) 4 + O(( ( 7) 6 ). (13) 

Depending on the parameters we may have either a 
second-order phase transition with b(T) > or a first- 
order phase transition with b(T) < 0. In the first case 
the transition temperature T$p is given by a(Tsp) = 
as 



A^l + e-^ 



2(E - £oo)e- 2 ^ + A s (l + e-" B -)e-^- 



(14) 



with (3 — > l/iksTsp). This transcendental equation 
takes a simple form in some limiting cases: 



E a 

Eoo 

A„ 



= 
E 



T, 



, r =A 2 g /(A g +E ) 



T SP = A g e Eo/ikBTsp) 
Tsp — > A„ — E + 2E 



(15) 



as a 



For illustration we present in Fig. |2| (a) Tsp 
function of the inter-chain coupling constant A g for 
E = 0.2meV. 
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FIG. 2. (a): Illustration of the spin-Peierls transition 
temperature Tsp for Eo = 0.2meV as a function of A 9 . 
(b): Phase diagram as a function of Eoo and Eo for fixed 
Tsp = 14.15X. Only the region Eo > Eoo is allowed (above 
the solid line). Above the dashed line the phase transition 
is of first-order, below it is of second-order. The crosses in- 
dicate the parameter values considered for comparison with 
CuGeO-i. 
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For b(Tsp) < we obtain a first-order phase transi- 
tion. This is the case for values of Eq larger than a 
certain critical value of the soliton energy E c , which is 
determined from b(Tsp) — as 



transition, as can be deducted from the closeness of the 
experimental points in Fig. |§| to the critical curves where 
E w E c (£00), compare Fig. |(b). 



E c — E n 




1 + 3 e JWT SP + 3e 2B 00 /T Sf 

6(l + e £ ~/^p) 2 



(16) 



For a fixed transition temperature T$p = 14.15K and 
using Eq. (|l4|), E c and the corresponding inter-chain 
coupling energy A g can be calculated as a function of 
Eqo. The resulting phase diagram is given in Fig. |^(b). 
The numerical results presented throughout this paper 
are obtained within the second-order regime, indicated 
by the crosses in Fig. ||(b). 



V. SELF-CONSISTENCY EQUATION 



The effective Hamiltonian Eq. contains three free 
parameters, namely -Boo, £-0 and A g . We examine two 
scenarios. The first is the case of a c > a with E^ = 0. 
The second is the case of a c < a « 0.35 with Eoo = 
0.15meV corresponding to a gap in the disordered phase 
of 2Boo = 0.3meVO. For each case we consider a range 
of Eq (see Fig. §(b)) within the second-order regime, 

Eoc < E < E c , 



Er. 



OmeV; E = 0, 0.1, 0.2, 0.277 meV 



(17) 



2E oc = 0.3 meV; E = 0.15, 0.3, 0.4, 0.474 meV (18) 

The largest value of E for each E^ corresponds to 
E^Eoo), compare Fig. ||(b). The experimental transition 
temperature T$p = 14.15K. of CuGeC>3 then determines 
the coupling constant AgEfl. 

The order parameter (cr) is determined self-consis- 
tently as a function of temperature by setting the deriva- 
tive of the free energy with respect to (cr) to zero, 



sinh(/?A g (cr)) 



' A„A 



-2/3 E s 



sinh 2 (/3A g (cr)) +■ 



-2f3E 3 



(19) 

In Fig. H we show the results for ~ (cr) 2 as a function of 
temperature for the parameters given in ( [l7|) and (18). 
We have also plotted the_measured intensity of an addi- 
tional super-lattice peako, which is proportional to the 
square of the lattice dimerization. We have normalized 
the experimental data such that agreement is obtained 
in the low-temperature regime. 

The comparison between theory and the data for 
CuGeC>3 does not lead to a determination of the frus- 
tration a but indicates closeness to a first-order phase 
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FIG. 3. The square of the spin-singlet order parame- 
ter ({cr) 2 , lines) as a function of temperature, for two dif- 
ferent values of E^, (soliton gap in the disordered phase) 
and several values of Eo (T — soliton gap). (a): 
E x = 0, i.e. E s = E (o) 2 . (b): 2Eoo = 0.3meV, i.e. 
E s = -Kepi + {Eq — Erx)(a) 2 . For comparison we plot the 
measurecO intensity of an additional super-lattice peak (hlled 
circles). 



VI. THERMODYNAMICS 

The entropy S and the specific heat cy are obtained 
from the free energy (Jfy via 



S 



dF 
df' 



T- 



df' 



(20) 
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In the disordered phase the entropy is (for E^ — 0) tem- 
perature independent with a value of fcsln(2)/2 per 
site, which is only half of the expected value for a spin 
1/2 chain. This is due to the fact that we did neglect up 
to now the spin-degrees of freedom of the domain-walls. 
In Fig. |J we present results for cy(T) for the parameters 
given by (|l7|) and (|l|). 
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FIG. 4. SpeciBc heat for different parameters as a function 
of temperature, (a): = 0. (b): 2E X = O.ZmeV. The 
inset in graph (a) shows the experimental dat£i of cy /T 
in units of [mJ/gK 2 ] versus T (10 -2 meV/K 2 corresponds to 
5.24mJ/gK 2 ). 

For small soliton excitation energies the results are typ- 
ically mean-field like. For values of Eo approaching the 
limit of the second-order phase regime the specific heat is 
strongly enhanced near T$p- It will diverge as the tran- 
sition becomes first-order. Note that the specific heat is 
linear in temperature for E^ < Eq < E c , in the limit 
T -^>Tsp and that the jump in the specific heat diverges 
as E -> E c like (E c - Eq)- 1 . Right at E = E c the 



specific heat diverges like (T$p — T) -1 / 2 . Note that a 
similar divergence ~ (T$p — T)~°' 4 |— has been reported 
in an early measurement for CuGeC>3c3 though the exact 
value of the specific-heat critical exponent for CuGeC>3 
is still controversialEZlEa 

In the inset of Fig. [| we present the measured mag- 
netic contribution to the specific heat of CuGeOsEi 
A straightforward comparison with the results of the 
mean-field theory is not possible since all the entropy 
of the effective Ising chain is released in a mean-field ap- 
proach. This corresponds to half of the entropy of the 
spin chain. Experimentalhz-jonly ~ 10% of the magnetic 
entropy is released at Tspc3, since the exchange constant 
J w 160K 3> Tsp and the measured specific heat is con- 
sequently smaller in magnitude than our mean-field re- 
sult. The neglect of the soliton dispersion relation (|^) 
is, on a microscopic level, the reason for this discrepancy 
between theory and experiment. A qualitative compar- 
ison is nevertheless possible and favors an Eq close to 
the first-order phase transition. 

Up to know we did not take the spin-degree of freedom 
of the solitons into account, as they just contribute a con- 
stant factor to the partition function in the paramagnetic 
case. As we have no magnetic interactions between the 
spins of different solitons in our model we can evaluate 
the magnetic susceptibility simply by Curie's law 



X(T) 



g 2 ^ B S(S+l) 
3k B T 



n(T) « 1.16-^1 n(T ), (21) 
k B T 



where fip = sW(2m e c) is the Bohr magneton, g = 2.15 
the measuredEj g- factor of the Cu 2+ — ion, S = 1/2 
and n(T) the density of thermally activated solitons 
per site. n(T) is obtained differentiating the free energy 
with respect to E s : 



i(T) = 



dF/L 



1 



-2/3 E s 



2A 



v /sinh 2 (^A g (a))+e- 



-2f3E 3 



(22) 



Above Tsp this expression reduces to 
n(T > T SP ) = ~ ' 



1 



(23) 



The mean number of solitons attains 1/4 per site, as 
the temperature goes to infinity, corresponding to half 
a soliton per dimer. The results for the magnetic sus- 
ceptibility are shown in Fig. ^. The fast drop of x(^) 
below Tsp for larger values of the soliton excitations en- 
ergies ispagain reminiscent of the experimental data for 
CuGeOaEj. Our susceptibility rises though much higher 
at Tsp than the experimental data which is a direct 
consequence of the neglected soliton dispersion (^). In 
the limit of large temperatures T 3> J the theoretical 
curve drops to 1/4 of the experimental points as a con- 
sequence of the aforementioned soliton density (23) (at 
300K « 2 J it has dropped to about 1/2 of the magnitude 
of the experimental data). 
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Note that for the curves plotted as solid lines in Fig. [|, 
Fig. |J, and Fig. [|, which are closest to the experimental 
data within the parameters chosen by us, entropy, specific 
heat, and susceptibilty are overestimated by about the 
same factor of five. 
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FIG. 5. Magnetic susceptibility for different parame- 
ters as a function of temperature, (a): — 0. (b): 
2Eoo — 0.3meV^ The inset in graph (a) shows the ex- 
perimental dataEj of x m units of [W~ 9 m 3 /mole] versus T 
(0Mfj,%/k B K correspond to 47.12 • 10" 9 m 3 /mole). 



VII. SINGLET-TRIPLET GAP 

The gap to triplet excitations is given by Eq. (^), 



A N = 2(E oo + {E -E oo ){v) 2 ) +c' J 



J 



2/3 



(24) 



with d = 2.33. A straightforward application of ( p4| ) 
would yield, compared with experiment, a much too large 



zero-temperature gap 2E + c'J(A g / J) 2 / 3 . This is so 
since the order parameter (ct)(T = 0) takes the value 
one in the molecular-field approximation, while it is much 
smaller than unity for CuGeC>3. We have therefore de- 
cided to use the parameter c! in (|24|) to fix A^(T = 0) 
to the experimentally observed value 2.5meV. The re- 
sults are giv£n in Fig. |^, together with the measured gap 
for CuGeOji 
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FIG. 6. The gap to triplet excitations for several pa- 
rameters as a function of temperature, (a): E^ — 0. (b): 
2Eoo = 0.3meV. 



VIII. CONCLUSIONS 

We have discussed a simple mean-field theory for spin- 
Peierls transitions applicable both for phonon-driven 
(a < ot c ) and for magnetically driven (a > a c ) spin- 
Peierls transitions. We have applied the approach to 
CuGeOa and found that it is not possible to determine 
uniquely from the experimentally measured temperature 



G 



dependence of the order-parameter the magnitude of the 
frustration parameter a. 

The theory allows both for a first-order and a second- 
order spin-Pcicrls transition, depending on the parame- 
ters of the model. We find that the parameters which 
fit experiments best indicate that CuGe03 is close to a 
first-order phase transition. 
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forward derivation of the S 2 ^ 3 scaling law. The support 
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